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Abstract We present a free energy functional for superfluid 3 He in the presence of 
homogeneously distributed impurity disorder which extends the Ginzburg-Landau 
free energy functional to all temperatures. We use the new free energy functional 
to calculate the thermodynamic potential, entropy, heat capacity and density of 
states for the B -phase of superfluid 3 He in homogeneous, isotropic aerogel. 

PACS numbers: 67.30.ef, 67.30.H-, 67.30.hm, 61.43.Hv 
Introduction 

Liquid 3 He infused into high porosity silica aerogel provides us with a model 
physical system for studying the effects of quenched disorder on an Fermi liquid 
with unconventional pairing ID. Indeed 3 He -aerogel provides us with a system 
where the effects of disorder can be explored over the entire range from weak to 
strong disorder without modifying the basic interactions of the Fermi liquid. 

Theoretical analysis by Thuneberg et al. 12 showed that the effects of scat- 
tering by a homogeneous, isotropic medium enhanced the stability of the Balian- 
Werthamer (BW) phase relative to the anisotropic ABM or Planar phase. These 
authors also calculated the reduction in the transition temperature, T c , the order 
parameter, 4(p), and condensation energy, AQ(p,T), in the Ginzburg-Landau 
(GL) limit. Experimentally, the equilibrium phase of 3 He-aerogel is consistent 
with the Balian-Werthamer (BW) state, modified by de-pairing, over the entire 
phase diagram, except perhaps in a narrow region of temperatures near T c at high 
pressures [3 4 5]. Measurements of the heat capacity of 3 He-aerogel at intermedi- 
ate and high pressures provide quantitative measurements of de-pairing, including 
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the reduction in AC/y H T CQ as well as indirect evidence of gapless excitations with 
a finite density of states, N(0), at the Fermi energy J5). 

We report results based on a free energy functional which extends the GL 
theory for superfluid 3 He in the presence of homogeneously distributed impu- 
rity disorder to all temperatures. This functional is obtained by a reduction of the 
Luttinger-Ward functional to leading order in the expansion parameters of Fermi- 
liquid theory. Results for the thermodynamic potential, entropy, heat capacity and 
density of states for the B-phase of superfluid 3 He-aerogel in the weak-coupling 
limit are reported. For impurity scattering in the unitary limit the B-phase exhibits 
gapless behavior over the entire phase diagram for high-porosity aerogels with an 
elastic mean- free path 4i = 1500 A. 

Free Energy Functional 

Serene and Rainer (SR) obtained a functional of the low-energy, quasiclassical 
propagator, g(p,£„), and self-energy, a(p,£„), for a superfluid Fermi liquid start- 
ing from the Luttinger-Ward-DeDominicis-Martin (LWDM) functional Q [8] of 
the exact propagator and self energy and subtracting the functional evaluated at 
the stationarity condition corresponding to the normal Fermi liquid [9|. This sub- 
traction confines the functional to the low-energy states near the Fermi surface^] 
Note that a, represents corrections to the leading order normal-state self-energy. 
Similarly, the subtracted ^-functional can be expressed as a functional of the low- 
energy propagator and effective interactions represented by block vertices that 
sum to all orders the bare interaction and higher order processes mediated by the 
high-energy propagator. So long as the block vertices do not introduce a new low- 
energy scale, then the contributions to the subtracted <J>-functional can be classi- 
fied by their order in one or more small parameters of Fermi-liquid theory. The 
SR functional has the form, 

AQ[g,8] = - l 2 S P '{ag} + An ln + <P [g] , (1) 
AQ la = ~\sp'{J d£, v $n® (-ie n % + % p + a)-ln®(-ie n % + Z p )}},(2) 

where Sp'{...} = 2N f Jd 3 R j ^ T^Te* Tr 4 {. . .} and Tr 4 {...} is the trace in 
Nambu space. The SR functional is general enough to describe inhomogeneous 
equilibrium states of a superfluid Fermi liquid which vary slowly on the atomic 
scale. In such cases the convolution product attached to the ln-functional is de- 
fined as a® fe(p,R) = expj| (V R • V* - Vp* • V R ) j a(p,R)%,R). For example 

<^ p ®a = (£ p — l\f • Vr)cT(p,R). The SR functional is stationary with respect to 
variations of g and a. The latter condition generates, 

g(p,R;£„) = f d$ p (;?,£„ , (3) 



1 We follow the notation in Ref. |9|, except that we denote 4x4 matrices in Nambu space by 
a widehat, e.g. g, while narrow hats refer to unit vectors, e.g. p, is the unit vector in the direction 
of the Fermi momentum, p/- = p/p. 
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which can be transformed a la Eilenberger into the quasiclassical transport equa- 
tion fflOl , i%e n — a, g +i\f-V R g = 0, and normalization condition, g(p,R;£„) 2 = 

— 7T 2 1. Stationarity of the SR functional with respect to variations of the propaga- 
tor, g — > g + 5 g, identifies the self energy with the skeleton expansion obtained 
from the "^-functional, 

8<P [g] = <P [g + 8g] -<P[g] = l 2 S P '{a5g} . (4) 

To leading order in the small expansion parameter the <£>-functional is given by 

& s = N f J d 3 R J d ®» d °*' v(p,pO T £f (p, R; e n ) ■ T £ f (p', R; , (5) 

£ " £ 'n 

where v(p • p') is the odd-parity, spin-triplet pairing interaction, while f (p, R; £„) 
and f(p, R; 8„) are the spin-triplet components of the the off-diagonal "anomalous" 
components of the propagator. This functional generates the mean-field pairing 
self energy, 

A(p,R)= J ^v(p,p')r£f(p',R;<). (6) 

The Matsubara sums include a cutoff, Q c , that restricts the pairing interaction to 
the low-energy band near the Fermi level. At this same order the diagonal mean- 
fields, i.e. the Landau molecular fields, also contribute. These terms are important 
for states perturbed by external fields or for inhomogeneous phases, but vanish 
for the homogeneous equilibrium state. For the applications discussed below the 
Landau molecular field self energy is omitted. 

In order to obtain a generalization of the GL functional we invert the self- 
consistency condition for the anomalous propagator to reduce the SR functional 
to a functional of the order parameter, A. Thus, <P n can be evaluated to give, 

<P a =N f ld 3 R^-j d ^\A(p 1 R)\ 2 , (7) 

for pure p-wave pairing. Similarly, in Eq. |IJ), — ^Sp' {Bg} = — 24> ff . 

In order to include the interaction of liquid 3 He with the structure of silica 
aerogel we explicitly include a one-body interaction, V = Vl, representing the in- 
teraction of quasiparticles with the static impurity potential, V(p,p'; {R,}). The 
impurity potential enters via the ln-functional, and the stationarity condition with 
respect to a generates Eq. d31), but now with the impurity potential as an addition 
to the self energy a. The impurity potential depends on the static configuration of 
the impurities, {R/}, which are treated as random variables governed by a proba- 
bility distribution, P({R,}). The impurity-averaged propagator is defined by aver- 
aging the propagator for a specific configuration of impurities over the probability 
distribution for a particular configuration, 



0(p,R;e„) = | d% p j\J d 3 RiP({Ri}) (ie„%- $ P -A -v({R,})) J . (8) 
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The impurity average generates an additional self energy term describing the scat- 
tering of quasiparticles and correlated pairs by the random potential. The deriva- 
tion of the transport equation carries through for the ensemble averaged propaga- 
tor, which now takes the form 

+ ivf -V R = O. (9) 

For uncorrelated impurities the impurity self energy, to leading order in 1 /zEt, 
is determined by mean density of impurities and the self-consistent t-matrix for 
impurity scattering, 

a, mp (p,R;e„) = n imp t(p,p,R;e„) (10) 
t(p,p',R;e„) = w(p,p')+AT/ / ^M(p,p")fl(p",R;e«)t(p",p',R;e„) • (11) 

For the results presented here we consider isotropic impurities with dominant scat- 
tering in the s-wave channel, i.e. u = uq\. In this case the t-matrix reduces to, 

t(R;e„) = u (l+Nf(g(p,R;e n ))^t(R;e n )^ , (12) 

where (g(p,R;£„)) p is the average of the propagator over the Fermi surface. For 
normal 3 He the propagator reduces to g N = — i'7rT 3 sgn(e„), and the normal-state t- 
matrix is parameterized by the s-wave scattering phase shift, So = tan -1 [%Nf uq), 
and reduces to t N (e„) = sm5oe~' Si;So ' 13 , where s £ = sgn(e„). In this minimal 

model for aerogel scattering, the mean density of impurities n, = \J &£, and the 
scattering rate, 1 /t, for normal-state quasiparticles are fixed by the mean free path, 
I = VfT, and scattering cross-section, a = AnXj do where kf = h/pf is the Fermi 

wavelength and do = sin<5g is the dimensionless cross-section [ 12 1. 

In order to average the free energy functional in Eq. <[TJ over the impurity 
ensemble we must average the ln-functional in Eq. Q. This is facilitated by the 
integral representation for the ln-functional, 

- Tr 4 {ln 8 (-iEn% + +4(p,R) + V({R,-})) } 

= Tr 4 |/T 3 £°^ (a^-^p-^R)-^,}))^ 1 ! . (13) 



For the subtracted functional we can now carry out the ensemble average, and 
interchange the ^ p integration with the auxiliary integration over the Matsubara 
energy to obtain the ensemble averaged free energy functional, 




H,e n -A- a imp , g 
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wher e fl N (A ) = — /7TT,sgn(A) and g(p,R;A) is obtained from the solution to the 



Eqs. ( 9 10 and 1 2 1 for fixed A (p, R) and auxiliary Matsubara energy, iX . This free 
energy functional can be used to evaluate the thermodynamic potential for a broad 
range of equilibrium states of unconventional superconductors in the presence of 
impurity disorder; it also represents an extension of the Eilenberger's free energy 
functional for conventional superconducting alloys in the Born impurity scattering 
limit [ 13 1. As an application we calculate the thermodynamic potential for the BW 
state in the presence of homogeneous, isotropic impurity scattering as a model for 
the B-phase of 3 He in silica aerogel. 

For spatially uniform states the solution to the transport equation and normal- 
ization condition for the propagator can be expressed as 

Q = -n—F^ =, (15) 



e2 + |A(p; £n )|2 

where e„ = £„ + in, mp t 3 (e„) and A = A(p) + A imp (e„) are the Matsubara energy 
and order parameter renormalized by the impurity self-energy. However, A imp (e„) 
vanishes for pure s-wave scattering because the order parameter vanishes when 
averaged over the Fermi surface, i.e. (A(p))« = III II . Thus, we obtain from the 
T 3 -component of the t-matrix, 



e„ = e n + -( ~ "; A " ' ' r 71 - : ■ (16) 




The self-consistency equation for the order parameter is obtained by evaluating 
Eq. ^ with the solution for f(p;£„) obtained from Eq. (15i. For unitary spin- 
triplet states (A x A* = 0) the gap equation becomes, 

J 47t t V^? + |A(p')l 2 

For pure 3 He the linearized gap equation (LGE) with A — > and 1/2t — > de- 
termines the bulk transition temperature, T CQ , in terms of the p-wave pairing inter- 
action, vi, and the bandwidth of low-energy states, Q c , i.e. K(r CQ ) = — 1/vi, with 
K(T) = 2nTY! £n > {l /£«) - ln(1.13X2 c /T). The LGE also regulates the Matsub- 
ara sum and removes the pairing interaction and cutoff in favor of the pure 3 He 
transition temperature, T CQ . 

In (f) A(p) = 2kT £ ( ^=== 2 f] A(p) . (18) 

\ T coJ „T VV e H + l A (P)| 2 £ "J 



The renormalization of the Matsubara energy encodes the effects of de-pairing 
by the scattering of pair-correlated quasiparticles off the distribution of impurities. 
This is evident in suppression of the superfiuid transition and the order parameter. 
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The suppression of the superconducting transition is given by the solution of Eq. 
( 18 1 for | A| -S- + . In this limit, e„ -)■ e„ + ^sgn(e„) and T — > T c , 



In 



'c 




(19) 



and we obtain the Abrikosov-Gorkov formula [2|. Note that the pair-breaking pa- 
rameter, a = 1 /2mT CQ = §a(p)/t, where = f /InkgT^ is the pure 3 He co- 
herence length and i = v/T is identified with the transport mean free path. The full 
range from weak to strong pair breaking can be explored by varying the pressure 
and aerogel density. The critical point in the phase diagram where T c = is given 
by a c ~ 0.281. For 98% aerogel with a mean free path of I ~ 1400 A this gives a 
critical pressure of p c a 6 bar, below which the superfiuid transition is completely 
suppressed |[T4l . For superfluidity to survive down to p = Obar the mfp must ex- 
ceed « 3000 A, while for £ < 650A superfluidity is completely suppressed for all 
pressures. 




Fig. 1 (Color online) Suppression of the B-phase order parameter in the unitary limit as a func- 
tion of a = \/2k%T cq . A(T) is normalized by Aq(Q) = \.T&k.BT C(j for pure 3 He-B in the weak- 
coupling limit (red curve). The curves for < a < 0.275 are shown for a in steps of 0.025. 



For the B-phase, |A(p)| 2 = \A\ 2 is isotropic, the angular integration over the 
Fermi surface is trivial and we can solve the gap equation numerically to obtain 
the magnitude of the order parameter. Figure [TJ shows the de-pairing effect on both 
the transition temperature and the magnitude of the order parameter over the full 
temperature range as a function of a for unitary scattering. 

The free energy functional for spatially uniform phases of 3 He-aerogel is ob- 
tained by carrying out the integration in Eq. ( |14| l over the impurity renormalized 
propagator in Eq. ( 15 I. For unitary spin-triplet states and impurity scattering in the 
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unitary scattering limit we obtain, 



Ati[A,A*] =N f v( |A(p)| 2 +27rr^ 



vi 



x 1 



1 

2r|e„| 



2t' 



In 



e„ - \/^ + |A(p)| 2 
Ve,? + |A(P)P N 



En 



i VWTJmW 



2t 



, (20) 



(21) 



This functional reduces to the weak-coupling functional for pure 3 He (Eq. 5.16 
of Ref. [9]) for T — > °°. One can also verify that the stationarity condition for the 
ensemble averaged functional, 8A£2/8A(ji)* = 0, generates the impurity renor- 
malized gap equation, i.e. Eq. ( fT7j ). 

The thermodynamic potential, entropy and heat capacity of BW phase with ho- 
mogeneous, isotropic disorder are obtained by evaluating the ensemble averaged 



free energy functional with the self-consistently determined solutions of Eqs. ( 18 1, 
d20l and d2T 




T/T co 

Fig. 2 (Color online) Thermodynamic potential (left panel) and entropy (right panel) for the B- 
phase of 3 He-aerogel in the unitary limit as a function of T/T c . for pair breaking ranging from 
< Of < 0.275 in steps of 0.025. Results for pure 3 He-B are shown as red curves. 



Thermodynamic Functions 

Figure [2] shows the thermodynamic potential (left panel) and entropy (right panel) 
for the dirty B-phase for values of the pair breaking parameter spanning the range 
from weak to strong pair breaking. Note that we plot the full thermodynamic po- 
tential, Q S (T) — Q N (T) + AQ(T), which includes the normal-state contribution 
X2 N (r) = — j/nT 2 , and we normalized the potential in units of ^y N T^ o , where 
y N = |7T Nf is the Sommerfeld coefficient for the normal state. The entropy was 
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Fig. 3 (Color online) Specific heat capacity normalized by fnT CQ (left panel) and by y N J (right 
panel) for the B-phase of 3 He-aerogel in the unitary limit as a function of T/T C(j for pair breaking 
ranging from < a < 0.275 in steps of 0.025. Results for pure 3 He-B are shown as red curves. 



computed by numerically differentiating the thermodynamic potential, S(T) = 
—d£2/dT\ w . For pure 3 He-B the entropy vanishes exponentially for T CT C|] , but 
for the disordered B-phase in the unitary limit the entropy vanishes linearly for 
sufficiently low temperatures, S(T) — > y s T, where y s is the Sommerfeld coeffi- 
cient for the gapless B-phase. The gapless regime is achieved for T <C %„ p , where 
Ymt — a/4 (0)/2t is the bandwidth of the impurity-induced excitation spectrum 
near the Fermi level. 



The heat capacity provides both a quantitative measure of the order that de- 
velops at T c , as well as the spectrum of gapless states near the Fermi level. Figure 
[3] (left panel) shows results for the specific heat, C v = TdS/dT\ v over the full 
temperature range. The reduction of the heat capacity jump at T c computed from 
the thermodynamic potential agrees with the value calculated from the GL theory 
of Ref. 121 to three significant figures for all values of the pair breaking param- 
eter. The right panel is the same data for the heat capacity plotted as C v /YnT vs. 
T /T CQ . For weak to modest pair-breaking, a = 1 /2mT CQ < 0.1, corresponding to 
4i > 1600A at p = 33bar, C v /T is a non-monotonic function of T, however the 
limiting low-temperature heat capacity is linear in temperature, C v — > YsT, where 
y s is the Sommerfeld coefficient for temperatures below the impurity bandwidth, 
T < %„ p . The non-monotonic behavior of C v /T for a < 0.1 results from the gap 
between the impurity band and the continuum states with £ > \A\ (see Fig. Q. 
The limiting value provides the Sommerfeld coefficient for the gapless B-phase, 
Ys = Yn^ m T^oC v /YtiT. One can infer from the calculated Sommerfeld coefficient 
for the gapless B-phase the density of states at the Fermi level, i.e. y s = I^A^O). 
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Density of States 

On the other hand, the density of states can be calculated directly by analytic 
continuation of the Matsubara propagator to real energies, 

3 (ie n ) ► /(e) = -n \~ R r w ' (22) 



? = H#4«^!?. (23) 

W(£)=« / Ullm/( £ )J, ,24) 

where -/V(e) is the density of states (DOS) at excitation energy £ measured relative 
to the Fermi level. Figure |4] shows the results for the DOS (left panel) from weak 
to strong pair breaking calculated in the unitary limit for T — > 0. Note that the 
excitation energy for each value of pairbreaking is normalized by the appropriate 
value of the order parameter, 4(0). The divergence of N(e) at £ = ±\A\ is sup- 
pressed by impurity scattering, and gapless excitations are created near the Fermi 
level for any a 7^ in the unitary scattering limit. 




Fig. 4 (Color online) Density of states as a function of excitation energy from weak to strong 
pair breaking (left). Comparison of the normalized density of states at the Fermi level (right 
panel), N(0)/Nf obtained by direct calculation from Eq. (24), and the Sommerfeld ratio, y s /yk 
computed from the limit, limr^oCs/r. 



The normalized density of states at £ = can be compared with that inferred 
from the calculated from the Sommerfeld coefficient. This provides a strong test 
of our calculations for the thermodynamic potential, entropy and heat capacity, 
as well as the validity of the ensemble averaged free energy functional obtained 
in Eq. (20 1. As Figure [4] shows the results for N(0)/Nf and y s /y N (right panel) 
are in perfect agreement. Note that N(0) /Nf and J s /y H vanish for a = and both 
approach 1 continuously for a — > (X c ~ 0.281. 
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In summary we have obtained an ensemble averaged free energy functional 
for the superfluid 3 He in in the presence of homogeneous, isotropic impurity dis- 
order. This functional is based on a quasiclassical reduction of the Luttinger-Ward 
functional. Results for the thermodynamic potential, entropy and heat capacity 
were obtained by numerical evaluation of the stationarity conditions. The pre- 
dicted Sommerfeld ratio based on the new functional agrees exactly with a direct 
calculation of the DOS from a solution of the quasiclassical transport equation. 
A quantitative comparison between theory and experimental measurements of the 
heat capacity of superfluid 3 He-B in aerogel [6| made at intermediate and high 
pressures will be discussed in a separate report. 
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